**Volume Title** 

ASP Conference Series, Vol. **Volume Number** 

**Author** 

© **Copyright Year** Astronomical Society of the Pacific 



Nature, formation and evolution of High Mass X-ray Binaries 

Sylvain Chaty 

O '. Laboratoire AIM (UMR 7158 CEA/DSM-CNRS-Universite Paris Diderot), 

CN ' Irfu/Service d'Astrophysique, CEA-Saclay, FR-91191 Gif-sur-Yvette Cedex, 

France, chaty@cea.fr 

Abstract. The aim of this review is to describe the nature, formation and evolution 

of the three kinds of high mass X-ray binary (HMXB) population: /. systems hosting 

Be stars (BeHMXBs), //. systems accreting the stellar wind of supergiant stars (sgH- 

Py ' MXBs), and ///. supergiant stars overflowing their Roche lobe. There are now many 

hM ' new observations, from the high-energy side (mainly from the INTEGRAL satellite), 

• , complemented by multi-wavelength observations (mainly in the optical, near and mid- 

i-G ■ infrared from ESO facilities), showing that a new population of supergiant HMXBs 

r^'^ has been recently revealed. New observations also suggest the existence of evolution- 

Q , ary links between Be and stellar wind accreting supergiant X-ray binaries. I describe 

$-H ' here the observational facts about the different categories of HMXBs, discuss the dif- 

tj^ _ ferent models of accretion in these sources (e.g. transitory accretion disc versus clumpy 

d . winds), show the evidences of a link between different kinds of HMXBs, and finally 

compare observations with population synthesis models. 



m 

p\J . 1. Introduction 

p 

[~^ . Nearly 50 years afte r the discovery of the first extra-solar X-ray source (ScoX-1; 

O I iGiacconi et alj [19621) . X-ray astronomy has now reached its age of reason, with a 

plethora of telescopes and satellites covering the whole electromagnetic spectrum, ob- 
taining precious observations of these powerful celestial objects. 

High energy binary systems are composed of a compact object - neutron star (NS) 
or black hol e (BH) - orbiting a companion star, and accreting matter from it (cf e.g. 



X 



;_] ■ IChatyl 12006 . for a review). The companion star is either a low mass star (typically 



P^ I 1 Mq or less, with a specti^al type later than B, called in the following LMXB for "Low- 

Mass X-ray Binary"), or a luminous early spectral type OB high mass companion star 
(typically > 10 Mq, called HMXB for "High-Mass X-ray Binary"). 300 high energy 
binary systems are known in our Galaxy: 187 L MXBs and 114 HMXBs (respectively 



62% and 38% of the total: iLiu et al]l2006L 120071) . 



Accretion of matter is different for both types of sourceo In the case of LMXBs, 
the small and low mass companion star fills and overflows its Roche lobe, therefore 
accretion of matter always occurs through the formation of an accretion disc. The com- 
pact object can be either a NS or a BH, Sco X-1 falling in the former category. In 



'I will not review here the so-called y-ray binaries, where most of the energy is in the GeV-TeV range, 
coming from interaction between relativistic electrons of the pulsar and the wind of the massive star. 
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the case of HMXBs, accretion can also occur through an accretion disc, for systems 
in which the companion star overflows its Roche lobe; however this is generally not 
the case, and there are two alternatives. The first one concerns stars with a circum- 
stellar disc, and here it is when the compact object - on a wide and eccentric orbit - 
crosses this disc, that accretion periodically occurs (case of HMXBs containing a main 
sequence early spectral type Be III/IV/V star, rapidly rotating, called in the following 
BeHMXBs). The second case is when the massive star ejects a slow and dense wind 
radially outflowing from the equator, and the compact object directly accretes the stellar 
wind through e.g. Bondy-Hoyle-Littleton processes (case of HMXBs containing a su- 
pergiant I/II star, called sgHMXBs). We point out that there also exists binary systems 
for which the companion star possesses an intermediate mass (typically between 1 and 
10 Mq, called IMXB for "Intermediate-Mass X-ray Binary"). 



2. Different types of High Mass X-ray Binaries 

2.1. Be X-ray binaries 

BeHMXBs usually host a NS on a wide and eccentric orbit around an early spectral type 
B0-B2e (with emission lines) donor star, surrounded by a circumstellar "decretion" disc 
of gas created by a low-velocity and high-density stellar wind of ~ 10^^ Mo/yr. This 
disc is characterized by an Ha emission line (whose width is correlated with the disc 
size) and a continuum fr ee-free or free-bound emission (causing an infrared excess: ICoa 
20001 : lNegueru elall2004ln . They exhibit transient and bright X-ray outbursts: /. "Type 



I" are regular and periodic each time the NS crosses the decretion disc at periastron: 
//. "Type IF' are giant outbursts at any phase, with a dramatic expansion of the disc, 
enshrouding the NS: Hi- "Missed" outbursts exhibit low He emission (due to a small 
disc or a centrifugal inhibition of accretion), iv. "Shifting outburst phases" are likely 
due to the rotation of density structures in the circumstellar disc. 

There are ~ 50 BeHMXBs in our Galaxy, and > 35 in the Small Magellanic 
Cloud (SMC). This large number in the SMC, instead of ~ 3 predicted by the galac- 
tic mass ratio 'smc"'* ~ ^^' ^^ li kely due to the bridg e of material between our 



Galaxy and the Magellanic Clouds (iMcBride et al.l 120081) . There is a large number 



of Supernova Remnants (SNRs) o f similar age (~ 5 M yr) suggesting an increased star- 



birth rate due to tidal interactions (IStanimirovic et al.l ll999). since the previous closest 



; (~ 3 M yr 
ic et al.| [l9 



SMC/Large Magellanic Cloud (LMC) approach was ~ 100 Myr ago, and new mas- 
sive stars probably formed the current HMXB population. A strong spatial correla- 
tion exists between emission line stars and stars aged 8 - 12 Myr, with BeHMXBs 
in the SMC (Mevsson nier & Azzopardilll993l : iMaragoudaki et al.ll200lh . The num- 
ber of HMX Bs is ther efore an indicator of star formation rate and starburst activity 
(iPopov et al.| [l998 : Gri mm et al.l 2003), and the Magellanic Clouds provide an uniform 
sample of BeHMXBs in a compact region, all at the same distance. 



^For more de t ails a bout BeHMXBs, I warmly recommend the excellent reference chapter by 
ICharles & Cod ( l2006h on optical/infrared emission of HMXBs. 



2.1.1. Formation of Be X-ray binaries 

BeHMXB populations are similar in th e Milky Way and SMC/LMC, when examinin g 



the spectral type of the companion star (Negueruel a & Coell2002l : lMcBride et al.ll2008h . 



However, while comparing the spectral type of Be stars hosted in X-ray binaries versus 
isolated Be, all located in our Galaxy, it appears that the spectral type distribution of 
Be in X-ray binaries is much narrower, beginning at 09 V (~ 22 Mq), peaking at BO 
V (~ 16 Mq) and stopping at B2 V (~ 8 - 10 Mq). The mass distribution of Be stars, 
ranging between 8 and 22 Mq, is consistent with the fact that wide orbits are vulnerable 
to disruption during SN event, especially for less massive B stars. The low mass range 
therefore supports the existence of SN kick velocities and angular momentu m loss, low 



mass stars being ejected, while heavier systems become sgHMXBs (Porte gies Zwart 
19951 : iMcBride et alJl2008h . 



The for mation of BeHMXBs has been exp lained by the so-called model of "re- 
juvenation" dRappaport & van den Heuvellll982h . by being a product of binary evolu- 
tion. In this model, the mass transfer speeds up first the rotation of the outer layers 
of the B spectral type secondary star, and then the star itself, which starts to turn so 
quickly that a circumstellar disc of gas is created, giving birth to the Be phenomenon 
(i.e. th ese stars are neither bom as fast ro tators, nor spun-up in the final main-sequence 
stages; ICod l2000l : ICharles & Cod 12006'). These systems are thus the result of mod- 



erately massive binaries undergoing a semi-conservative mass transfer evolution, with 
wide orbits (200 to 600 days) produced before the first SN ev ent of the system, and 



eccentric orbits due to s mall asymmetries during the SN event ('van den He u veil [19831 : 



Verbunt & van den Heuvel 1995) . This model is consistent with the spectral types of 



the companion stars between 09 and B2V: stars less massive than ~ 8 Mq are ejected 
during the SN event, and stars more massive than ~ 22 Mq become supergiant stars. 



2.1.2. Structure of Be circumstellar disc 

In BeHMXBs, Be circumstellar discs are naturally truncated due to the presence of the 
NS, as a result of tidal torques at certain resonance points (where the Keplerian period 
Pk is an integer fraction of the orbital period Porb)> thus preventing any accretion of 
matter beyond these points. The accretion on the NS is unlikely for BeHMXBs where 
the NS is on a circular orbit, with a disc always truncated at a fixed size, smaller than 
the Roche lobe. This creates a persistent low-level X-ray emission from the stellar 
wind, and occasional Type II outbursts. Instead, BeHMXBs with a NS on a highly 
eccentric orbit will allow for periodic accretion, with the size of the rotating truncated 
disc depending on the orbital phase, and at periastron the disc can include the N S orbit, 
thus creating Type I X-ray outbursts (lOkazaki & Negueruelall200]1 :lReig"2007'). 



In addition, recent observations show that circumstellar discs in BeHMXBs also 
exhibit cycles of activity, likely due to their formation and dispersion. BeHMXBs there- 
fore exhibit 3 periods: a spin period Pspin> an orbital period Porb and a super-orbital 
period Psup- The first BeHMXB to have s hown this cycle of activity is A 0538-66, with 
Porb - 16.65 days and Pjup - 421 days (iMcGowan & CharlesI 120031) . Then, a com- 



prehensive study of the MACHO and OGLE light curves over 1 8 years light-curves of 
~ 20 BeHMXBs show that Psup, extending fr om 500 to 4000 days, seems to be corre- 
lated to Porb, from 10 to 300 days (Rajoelima nana et al.ll201lh . While the origin of this 



variation is still unknown, this correlation shows that it is certainly due to binarity, and 
not intrinsic to the Be star. 
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2.2. Supergiant X-ray binaries 

These systems usually host a NS on a circular orbit around an early spectral type su- 
pergiant OB donor star, with a steady wind outflow. They are separated in two distinct 
groups: Roche lobe overflow and wind-fed systems. The former group constitutes the 
classical bright sgHMXBs with accreted matter flowing via inner Lagrangian point to 
the accretion disc, causing a high X-ray luminosity (Lx ~ 10^^ ergs^^) during out- 
bursto The later group concerns close systems (forb < 15 d) with a low eccentric- 
ity, the NS accreting from deep inside the strong steady radiative and highly super- 
sonic stellar wind. This creates a persistent X-ray emission at regular low-level effect 
{Lx ~ 10^^~^^ erg s'^), exhibiting rare Type II outbursts, and no Type I outbursts. These 
systems exhibit large variations on short timescales, due to wind inhomogeneities. Dur- 
ing their long term evolution, the orbits of sgHMXBs will tend to circularize mor e 
rapidly with time, while the rate of mass transfer steadily increases JKaper et aLll2004) . 



A milestone in the evolution of close binary systems takes place during what is 
called the "common envelope phase". This phase is initiated when the compact object 
enters inside the envelope of the companion star, in an orbit which is rapidly decreas- 
ing due to a large loss of orbital angular momentum. This phase has been invoked by 



PaczynskJI (Il976h to explain how high energy binary systems with very short Porb can be 



formed, while both components of these systems - large stars at their formation - would 
not have been able to fit inside a binary system with such a small orbital separation. This 
phase of inward spiralling, currently taken into account in population synthesis mod- 
els, but never observed yet, probably because it is short (m odels predict a maximum 
duration of common envelope phase of only ~ 1000 years; iMeurs & van den Heuvel 



I989I) compared to the lifetime of a massive star (~ 10^ ' years), is an ingredient 



of prime im portance for understanding the evolution of high energy binary systems 



dTauris & van den Heuvel 2006) . 



2.3. The Corbet diagram 

In case of an accreting NS in an HMXB, it is possible to detect a pulsation, correspond- 
ing to the NS spin period /"spin- These X-ray accretion-powered pulsars are divided ac- 
cording to their nature and the dominant accretion process - BeHMXBs (disc-accreting 
systems) or sgHMXBs (wind-accreting and Roche lobe overflow systems) - in dis- 
tinct locations in the Corbet diagram, representing the NS Pspin versus the system Porb 
(ICorbetlll98 d. see Figure[T]l. This diagram is a valuable tool to study the interaction and 



feedback between the NS and accreted matter, and the influence of the local absorbing 
matter, the location of the different systems being determined by the equilibrium period 
reached by the rotation of the NS accreting matter on its surface. 

It is clear from this diag ram that a strong correlation exists for BeHMXBs - 
^spin "^ (^orb)^ (ICorbellll984l) -, due to the efficient transfer of angular momentum 



when the NS accretes material from the Be star decretion disc. Intuitively, a small 
(wide) orbit presents on average a high (low) stellar wind density, increasing (decreas- 
ing) the accretion pressure, and thus accelerating (slowing) the NS rotation, which in 
turn increases (decreases) the centrifugal inhibition, preventing (allowing) more accre- 



dit is interesting to note that CygX-1 is the only sgHMXB with Roche lobe overflow and stellar wind 
accretion hosting a confirmed BH. 



tion of matter. In contrast, the lack of correlation for sgHMXBs suggests that in this 
case, wind accretion is very inefficient to transfer angular momentum. 

High energy properties, and in particular the accretion efficiency, of a compact 
object accreting from a circumstellar disc, a dust cocoon or the stellar wind of the com- 
panion star, will strongly differ, depending on the geometry of the system, altogether 
with the size and shape of the optically thin and fully ionized region of gas, similar to 



3 



the Stromgren sphere (see e.g. |Pringlg|l973l : iHatchett & McCraylll977h . In particular. 



Pspin in HMXBs is determined by the stellar wind characteristics. sgHMXBs exhibit 
a spherically-symmetric radiation-driven wind, with a density p ex r~^ and a velocity 
600 - 900 km s~^ while BeHMXBs present an equatorial stellar wind with a density 



pocr ^ ^^ ^-^ dropping faster, and a lower velocity ~ 200 - 300 km s ^ (IWaters et all 



19881) . Therefore, larger density and velocity gradients at the distance of the NS make 
the wind-fe d accretion transfer of angular momentum more efficient in BeHMXBs than 
SgHMXBs (IWaters & van Kerkwiiklfl989h . 



Accretion can occur on a magnetized NS only if the pressure of the infalling mate- 
rial is greater than the centrifugal inhibition (corresponding to the Alfven radius located 
inside the magnetospheric boundary). The NS should then reach an equilibrium rotation 
period Peq corresponding to the corotation velocity Vc (at the magnetospheric radius) 
equal to the Keplerian speed Vk- If Vc > Vk (coiTcsponding to Pspin < Peq) the Pro- 
peller mechanism increases P spin by ejecting material, taking away angular momentum 
(llUarionov & Sunyaevlll975l) . If Vc < Vk (corresponding to /"spin > ^eq) the accretion 
eiffier reduces or increases Pgpin (spinning up or down), depending on the direction of 
angular momentum with respect to the NS spin. 

Taking into account the density of the stellar wind, and assuming a steady accretion 
rate with the angular momentum of same direction than the NS spin, /"spin should reach 



^eq °^ p~^^^ dWaters et al.lll988b . However, this is not the case, neither for BeHMXBs 
nor for sgHMXBs. Pgpin of BeHMXBs does not correspond to /"gq. constantly adj usting 
to the changing conditions in the wind, reflecting values of earlier evolution stage (King 
[1991). CuiTcnt Pspin of NS in sgHMXBs is longer than predicted, closer to Peq of a NS 
embedded in the st ellar wind, while the compani on star was still a main sequence star 



of spectral type O (IWaters & van Kerkwijkll 19891) 



3. The y-ray sky seen by the INTEGRAL observatory 

The INTEGRAL observatory is an ESA satellite launched on 17 October 2002 by a 
PROTON rocket on an eccentric orbit. It hosts 4 instruments: 2 y-ray coded-mask 
telescopes - the imager IBIS and the spectro-imager SPI, observing in the range 10 
keV-10 MeV, with a resolution of 12' and a field-of-view of 19° -, a coded-mask tele- 
scope JEM-X (3-100 keV), and an optical telescope (OMC). The y-ray sky seen by 
INTEGRAL is very rich, since 723 sources have been detected by INTEGRAL, reported 
in the A''^ IBIS/ISGRI so ft y-ray catalogu e, spanning nearly 7 years of observations in 



the 17-100 keV domain (iBird et al.ll2010ll I Among these sources, there are 185 X-ray 



binaries (representing 26% of the whole sample of sources detected by INTEGRAL, 
called "IGRs" in the following), 255 Active Galactic Nuclei (35%), 35 Cataclysmic 
Variables (5%), and ~ 30 sources of other type (4%): 15 SNRs, 4 Globular Clusters, 3 



'^See an up-to-date list at\http://irfu.cea.fr/Sap/IGR-Sources/\ maintained by J. Rodriguez and A. Bodaghee 
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Figure 1 . Corbet diagram (NS Pspin versus system Poib), showing the three HMXB 
populations: /. systems hosting Be stars (BeHMXBs), //. supergiant stars overflow- 
ing their Roche lobe, and ///. systems accreting the stellar wind of supergiant stars 
(sgHMXBs). The distinct location of these three types of HMXBs results from the 
interaction and equilibrium reached between accretion of matter and NS spin, al- 
together with the general characteristics of binary systems (such as mass ratio and 
orbital separation). The "misplaced" SFXT IGR J18483-0311 is indicated (credit 
J. A. Zurita Heras). 



Soft y-ray Repeaters, 2 y-ray bursts, etc. 215 objects still remain unidentified (30%). 
X-ray binaries are separated into 95 LMXBs and 90 HMXBs, each category repre- 
senting ~ 13% of IGRs. Among identified HMXBs, there are 24 BeHMXBs and 19 
SgHMXBs, representing respectively 31% and 24% of HMXBs. 

It is interesting to follow the evolution of the ratio between BeHMXBs and sgH- 
MXBs. During the pre-INTEGRAL era, HMXBs were rn ostly BeHMXBs. For in- 
stance, in the catalogue of 130 HMXBs bv lLiu et al.l fcOOO'), there were 54 BeHMXBs 
and 5 sgHMXBs (respectively 42% and 4% of the total number of HMXBs). Then, 
the situation changed drastically with the first HMXBs identified b y INTEGRAL: i n 



the catalogue of 114 HMXBs (-^128 in the Magellanic Clouds) of lUuetalJ (l2006h . 



60% of the total number of HMXBs were firmly identified as BeHMXBs and 32% as 
SgHMXBs. Therefore, while the ratio of BeHMXBs/HMXBs increased by a factor of 
1.5, the sgHMXBs/HMXBs ratio increased by a larger factor of 8. The ISGRI energy 
range (> 20 keV), immune to the absorption that prevented the discovery of intrinsically 
absorbed sources by earlier soft X-ray telescopes, allowed us to go from a study of indi- 
vidual SgHMXBs (such as GX 301-2, 4U 1700-377, VelaX-1, etc.) to a comprehensive 
study of the characteristics of a whole population of HMXBs... 

The most important result of INTEGRAL to date is the discovery of many new high 
energy sources - concentrated in the Galactic plane, mainly towards tangential direc- 
tions of Galactic arms, rich in star forming regio ns -, ex hibiting common characteristics 
which previously had rarely been seen (see e.g. IChaty & Filliatre 2005). Many of them 
are HMXBs hosting a NS orbiting an OB companion, in most cases a supergiant star. 
Nearly all the INTEGRAL HMXBs for which both Pspin and Porb have been measured 



are located in the upper part of the Corbet diagram (see Figure [T). They are X-ray 
wind-accreting pulsars typical of sgHMXBs, with longer pulsation periods and higher 
absor ption (by a factor ~ 4) compared to previously known sgHMXBs (Bodaghee et all 
l2007b . They divide into two classes: some are very obscured, exhibiting a huge intrin- 
sic and local extinction, - the most extreme example being the highly absorbed source 
IGR J163 18-4848 (Filhatre &Chatv 2004 ) -, and the others are sgHMXBs exhibit- 
ing fast and transient outbursts - an unusual characteristic a mong HMXBs -. Thes e 
are therefore called Supergiant Fast X-ray Transie nts (SFXT s,'Negu eruela et alj|2006h . 
with IGR J17544-2619 being their archetype (Pellizza et al.:.2006.) . 



3.1. Obscured HMXBs: IGR J16318-4848, an extreme case 

IGR J163 18-4848 was the first source di scovered by IBIS/ISGRI on INTEGRAL on 
29 January 2003 (ICourvoisier et alj l2003h . with a 2' uncertainty. XMM-Newton ob- 
servations revealed a comptonised spectru m exhibiti ng an unusually high level of ab- 
sorption: A^H ~ 1.84 X 102^_cnrfjMatt& Guain azzill2003l ). The accurate localisation 
by XMM-Newton allowed iFiUiatre & Chatyl (I2004r) to rapidly trigger ToO photomet- 
ric and spectroscopic o bservations in optic al/infrared, leading to the confirmation of 
the optical counterpart dWalter et all l2003b and to the discovery of the infrared one 
dFilliatre & Chatvll2004l) . The extremely bright infrared source (B > 25.4±1; I = 16.05± 
0.54; J ^ 10.33 + 0.14; H ^ 8.33 ± 0.10 and K^ = 7.20 ± 0.05 magnitudes) exhibits an 
unusually strong intrinsic absorption in the optical (Ay = 17.4 magnitudes), 100 times 
stronger than the interstellar absorption along the line of sight (Ay = 11.4 magnitudes), 
but still 100 times lower than the absorption in X-rays. This led iFilUatre & Chaty 
([2004) to suggest that the material absorbing in X-rays was concentrated around the 
compact object, while the one absorbing in optical/infrared was enshrouding the whole 
system. 

Near-infrared spectroscopy in the 0.95 - 2.5 jum domain allowed iFilliatre & Chaty 



(|2004|) to identify the nature of the companion star, by revealing an unusual spec- 
trum, with many strong emission lines: H and He I (P-Cyg) lines characteristic of 
dense/ionised wind at v = 400km/s, He II lines signatures of a highly excited region, 
[Fe II] lines reminiscent of shock heated matter, Fe II lines emanating from media of 
densities > 10^ - 10^ cm~^, and Na I lines coming from cold/dense regions. All these 
lines originate from a highly complex, stratified circumstellar environment of various 
densities and temperatures, suggesting the presence of an envelope and strong stellar 
outflow/wind responsible for the absorption. Only lum inous early spectr al type stars 
such as sgB[e] show such extreme environments, and IFilliatre & Chatyl (12004) con- 
cluded that IGR J 163 18-4848 was an unusual HMXB hosting a sgB[e] with charac- 
teristic luminosity 10^ L©, mass 30 M©, radius 2 Rn and tempera ture T - 20250 K, 



located at a distance between 1 and 6 kpc (see also lChaty & Filliatre.i2005l). This sour ce 
is therefore the second HMXB hosting a sgB[e] star, after CI Cam dClark et al.ll 19991) . 
By combining optical to mid-infrare d (MIR) observation s, and fitting these obser- 
vations withasimple spherical blackbody. JRahoui et al.ld2008h showed that IGR J 163 18- 
4848 exhibited a MIR excess, interpreted as being due to the strong stellar outflow em- 
anating from the sgB [e] companion star. They found that this star had a temperature of 
Ti, = 22000 K and radius 7?^^ = 20.4 Rq = O.l a.u., consistent with a supergiant star, 
with an extra component of temperature T = 1 100 K and radius R = 1 1.9 /?* = 1 a.u., 
with Ay = 17.6 magnitudes. Recent MIR spectroscopic observations with VISIR at the 
VLT and Spitzer showed that the source was exhibiting strong emission lines of H, He, 
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Figure 2. Broadband near- to mid-infrared ESO/NTT-H VISIR and Spitzer red- 
dened spectrum of IGR J 163 18 -4848, from 0.9 to 35 fim (from IChaty & Rahouil 

unni). 



Ne, PAH, S i, proving that the extra absorbing component was made of dust and cold gas 
(IChaty & R ahoui 201 1). By fitting the optical to MIR spectra with a more sophisticated 
aspheric disc model developped for HAEBE objets, and adapted to sgB[e] stars, they 
showed that the supergiant star was surrounded by a hot rim of dust at 5500 K, with a 
warm dust component at 900 K around it (Chaty & Rahoui 2011, see Figure O. 

By assuming a typical Porb of 10 days and a mass of the companion star of 20 Mq, 
we obtain an orbital separation of 50 /?©, smaller than the extension of the extra com- 
ponent of dust/gas {= 240 Rq), suggesting that this dense and absorbing circumstellar 
material envelope enshrouds the whole binary system, like a cocoon of dust (see Figure 
m left panel). We point out that this source exhibits such extreme characteristics that it 
might not be fully representative of the other obscured sources. 



3.2. Supergiant Fast X-ray Transients 

SFXTs constitute a new class of ~ 12 sources identified among the recently discovered 
IGRs. They are HMXBs hosting NS orbiting sgOB companion stars, exhibiting pecu- 
liar characteristics compared to "classical" HMXBs: rapid outbursts lasting only for 
hours, faint quiescent emission, and high energy spectra requiring a BH or NS accretor. 
The flares rise in tens of minutes, last for ~ 1 hour, their frequency is ~ 7 days, and 
their luminosity Lx reaches ~ 10^^ erg s~' at the outburst peak. SFXTs can be divided 
in two groups, according to the duration and frequency of their outbursts, and their 
yssi ratio. "Classical" SFXTs exhibit a very low quiescence Lx and a high variability, 

'^min 

while "intermediate" SFXTs exhibit a higher average Lx, a lower yssi and a smaller 
variability factor, with longer flares. SFXTs might appear like persistent sgHMXBs 



with < Lx > below the can onical value of ~ 10^^ ergs ', and superimposed flares 



(IWalte r & Zur ita Herasl l2007'). but there might be some observational bias in these gen- 
eral characteristics, therefore the distinction between SFXTs and sgHMXBs is not yet 
well defined. While the typical hard X-ray variability factor (the ratio between outburst 
flux and deep quiescence) is less than 20 in classical/absorbed systems, it is higher than 
100, reaching lO'* for some SFXTs (some sources can exhibit flares in a few min utes. 



like for instance XTEJ1739-302 & IGR J 17544-26 19, see e.g. iPellizza et al "]l2006h. 



To explain the emission of SFXTs in the context of sgHMXBs, iNegueruela et al 



(l2008h invoked the existence of two zones in the stellar wind created by the super 



giant star, of high and low clump density. This would naturally explain both the X-ray 
lightcurves, each outburst being due to the accretion of single clumps, and the smooth 
transition between sgHMXBs and SFXTs, and the existence of intermediate systems; 
the main difference between classical sgHMXBs and SFXTs being in this scenario 
the NS orbital radius. Indeed, a basic model of porous wind with macro-clumping 



(INegueruela et all 120081) predicts a substantial change in the properties of the wind 
"seen by the NS" at a distance r ~ 2 /?+, where we stop seeing persistent X-ray sources. 
There are 2-regimes: at r < 2 /?^^ the NS sees a large number of clumps, embedded in 
a quasi-continuous wind; and at r > 2 7?^ the density of clumps is so low that the NS 
is effectively orbiting in an empty space. NS in classical sgHMXBs can only lie within 
2 Ri, of the companion star. 

However there exists other models to explain these accretion processes, like the 



formation of transient accretion discs ( Ruffert 1996; Ducci et al.ll201Q) and the accre- 
tion with centrifugal/magnetic barriers (IBozzo et al.ll2008h . 

3.2.1. The intermediate SFXT IGR J18483-0311 

X-ray properties of this system were suggesting an SFXT nature (ISguera et al.ll2007r) . 



exhibiting however an unusual behaviour: its outbursts last for a few days (to compare 
to hours for classical SFXTs), and the ratio L^ax/Lmtn only reaches ~ 10^ (compared 
to ~ 10"* for classical SFXTs), therefore suggesting a high level quiescence. Moreover, 
^orb = 18.5 days is low compared to classical SFXTs with large and eccentric orbits. 
Finally, Porb and Pspin - 21.05 s values locate it inbetween Be and sgHMXBs in the 
Corbet diagram (see Figure [T]). Rahoui & Chatyl (i2008.) identified the companion star 
of this system as a B0.5 la supergiant, unambiguously showing that this system is an 
SFXT. Furthermore, they suggested that this system could be the first firmly identified 
intermediate SFXT, characterised by short, eccentric orbit (with an eccentricity e be- 
tween 0.4 and 0.6), and long outbursts: an intermediate SFXT nature would explain the 
unusual characteristics of this source among classical SFXTs. 

3.2.2. What is the origin of "misplaced" sgHMXBs? 



Liu et al.l (1201 ID noted that there are two "m isplaced" SFXTs in the Corbet diagram: 



IGR Jl 1215-5952 (a NS orbiting a Bl la star; Ne gueruela et al.ll2005n and IG R Jl 8483 



0311 (described in the previous paragraph; see Figure [T])- According to iLiu et al. 
(l201lh . these 2 SFXTs can not have evolved from normal main sequence OB -type 
stars, since /. their NS are not spinning at the equilibrium spin period Pgq of O V stars, 
//. they have not been able to spin up after reach ing Pga, and ///. th eir NS have not 
yet reached Peq of sgHMXBs (see e.g. IWaters & van Kerkwijki ll989l and section [23] ). 



They must therefore be the descendants of BeHMXBs (i.e. hosting 0-type emission 
line stars), after the NS has reached Pgq, suggesting that some HMXBs exhibit two pe- 
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riods of accretion, the first one as BeHMXB and tlie second one as sgHMXB (iLiu et al. 



20111) . This important result suggests that there must be many more such intermediate 



SFXTs. 

3.3. A scenario towards the Grand Unification of sgHMXBs 

In view of the results described above, there seems to exist a continuous trend, with the 
observed division between classical sgHMBs and SFXTs being natu rally explained b 



simpl e geometrical differences in the orbital configurations (see e.g. iChaty et al.ll200 
20Toh : 



Classical (or obscured/absorbed) sgHMXBs: These systems (like IGR J16318-4848) 
host a NS on a short and circular orbit at a few stellar radii only from the star, 
inside the zone of stellar wind high clump density (/?orb ~ 2 /?*). It is constantly 
orbiting inside a cocoon of dust and/or cold gas, probably created by the compan- 
ion star itself, therefore creating a persistent and luminous X-ray emission. The 
cocoon, with an extension of ~ 10 /?* = 1 a.u., is enshrouding the whole binary 
system (see Figure [3j left panel). 

Intermediate SFXT systems: In these systems (such as IGRJ18483-0311, Porb - 
18.5 days), the NS orbits on a short and circular/eccentric orbit outside the zone 
of high clump density, and it is only when the NS penetrates inside the narrow 
transition zone between high and low clump density, close to the supergiant star, 
that accretion takes place, and that X-ray emission arises, with possible periodic 
outbursts. 

Classical SFXTs: In these systems (such as XTEJ1739-302, Porb ^ 50 days), the NS 
orbits outside the high density zone, on a large and eccentric orbit around the 
supergiant star, and exhibits some recurrent and short transient X-ray flares, when 
it comes close to the star, and accretes from clumps of matter coming from the 
wind of the supergiant. Because it is passing through more diluted medium, the 
j^ ratio is higher and the quiescence lasts for longer time for classical SFXTs 
compared to intermediate SFXTs (see Figure [3l right panel). 

Although this scenario seems to describe quite well the characteristics currently 
seen in sgHMXBs, we still need to identify the nature of many more sgHMXBs to 
confirm it, and in particular Porb and the dependance of the column density with the 
orbital phase of the binary system. 



4. Formation and evolution of X-ray Binaries 

4.1. Galactic distribution of X-ray binaries 

The evidence that the stellar birthplace is very important in the formation and evo- 
lution of X-ray binaries was initi ally noticed ex amining the Galactic d istribution of 
X-ray sources detected by Ginga (iKoyama et al.|[l99Q) and RXTE/ASM dGrimm et al. 



20021) . Thereafter, a statistical analysis of a more compl e te sample of sources prov ided 



by INTEGRAL has bee n undertaken by lLutovinov et all dloO^ . lDean et all (l2005l) and 



Bodaghee et al.l (120071) . These works show that LMXBs are concentrated towards the 
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Figure 3. Scenario illustrating two possible configurations of INTEGRAL sources: 
a NS orbiting a supergiant star on a circular orbit (left image); and on an eccen- 
tric orbit (right image), accreting from the clumpy stellar wind of the supergiant. 
The accretion of matter is persistent in the case of the obscured sources, as in the 
left image, where the compact object orbits inside the cocoon of dust enshroud- 
ing the whole system. On the other hand, the accretion is intermittent in the case 
of SFXTs, which might correspond to a compact object on an eccentric orbit, as 
in the right im age. A 3D animation of these sources is available on the website: 
^ttp://www.aim.univ-paris7.fr/CHATY/Research/hidden.html 



central regions of the Galaxy, and more precisely in the Galactic bulge, and then grad- 
ually decrease outside, while HMXBs, underabundant in the central kpc, are spread 
over the whole Galactic plane, confined in the disc, exhibiting an uneven distribution, 
preferentially towards the tangential directions of the spiral arms. This spatial distri- 
bution was expected, because LMXBs are constituted of companion stars belonging to 
an old stellar population, located in the Galactic bulge where they had time to migrate 
out of the Galactic plane {\b\ > 3 - 5°). On the contrary, HMXBs - containing young 
companion stars - remain close to their stellar birthsite. 

To perform this study, one has to derive the distribution of Galactic HMXBs from 
the catalogues of sources detected by current high-energy satellites such as XMM, 
Chandra, Swift and INTEGRAL, taking into account the observational bias, some ar- 
eas being more observed than others. However, the main difficulty is that most of 
the distances of HMXBs are not accurately determined. To overcome this difficulty, 
a possibility is to determine the distance of all known HMXBs, by fitting the optical 
and infrared magnitudes with a black body corresponding to the companion star spec- 
tral type, and comparing their distribution with catalogs of massive stars, active OB 
stars and star form ation complexes (SFCs:, Russell. 2003.) . A preliminary result of this 
method, detailed in lColeiro & Chatyl (1201 ih . gives a typical clustering size of HMXBs 
with SFCs of 0.3 kpc, with an average distance between these clusters of 1.7 kpc, there- 
fore suggesting a correlation (see Figure |4|. 

The appearance of HMXBs should theoretically follow the pa ssage of the wave 
density associated with the rotation of the spiral arms dLin et al.lll969 ). inducing a burst 
of stellar formation. In this context, the models predict a time lag between star for- 
mation and HMXB apparition, depending on the position in the Galaxy, an d reflecting 



the ra nge of initial masses of both stars composing each binary system (iDean et al. 



20051) . However, since the formation of an HMXB only takes tens of million years 
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Figure 4. Galactic distribution of HMXBs (blue stars) - with distance derived by 
adjusting the spectral type on optical/infrared magnitudes of the companion star -, 
and of SFCs (green circles, symbol size proportional to SF C activity; Russeil 2003). 
The Sun is indicated by the red dot at (8.5, 0) (from lColeiro & Chatyll201 Ih . 



(iTauris & van den Heuvelll2006h . it is a good marker for the transition from density 
wave in a region of the Galaxy. For such a statistical study to be accurate, we not only 
need many HMXBs for which we accurately know the distance, but also a precise and 
reliable kinematic model of the spiral arm structure of our Galaxy (Coleiro & Chaty 
subm.). 



4.2. Link with population synthesis models 

sgHMXBs revealed by INTEGRAL will allow us to better understand and constrain 
the formation and evolution of X-ray binary systems, by comparing them to numerical 
study of LMXB/HMXB population synthesis models. For instance, these new sys- 
tems might represent a precursor stage of what is known as the "common envelope 
phase" in the evolution of LMXB/HMXB s, when the orbit has shrunk so much that 
the NS begins to orbit inside the envelope of the supergiant star. In addition, many pa- 
rameters do influence the various evolutions, from one system to another: differences 
in mass, size, orbital period, age, rotation, magnetic field, accretion type, stellar end- 
points, etc. Moreover, stellar and circumstellar properties also influence the evolution 
of high-energy binary systems, made of two massive components likely born in rich 
star forming regions. 

The formation and evolution of the supergian t INT EGRAL sources is mentioned in 
the very nice review about relativistic binaries by Ivan den Heuvel (2009). Porh in later 
evolutionary phases being linearly dependent of initial Porb, we can derive that cur- 
rently wide Porb ~ 100 days OB-supergiant Wr^G/ML binaries require initial systems 
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with Porb ~ 10 days dvan den Heuvelll2009r) . Systems having long forb are expected to 



survive the common envelope phase, and may then either end as close eccentric binary 
radio pulsar systems (double NS), or in some cases as BH-NS binaries. 

No such system is known yet, however some of these systems might harbour a BH 
as the compact object. Of course NSs are easier to detect through X-ray pulsations, but, 
as Carl Sagan pointed it out, "absence of evidence is not evidence of absence" . We 
should look for BHs orbiting supergiant companion stars in wind-accreting HMXBs, 
however this is only feasible through observational methods involving detection of ex- 
tremely faint radial velocity displacement due to the high mass of the companion star, 
or through extremely accurate radio measurements that will be available in the future. 
On the other hand, massive stars los e so much ma tter during their evolution that they 



might always finish as NSs (see e.g. iMaeder & Meynet.2008.) . If this is the case, then 



such systems hosting BHs might not form at all. 

Finally, these sources are also useful to look for massive stellar progenitors, for 
instance giving birth to coalescence of compact objects, through NS/NS or NS/BH 
collisions. They would then become prime candidates for gravitational wave emitters, 
or even short-hard y-ray bursts. 



5. Conclusions and perspectives 

Let us first recall the INTEGRAL legacy: 

• The INTEGRAL satellite has quadrupled the total number of known Galactic sgH- 
MXBs, constituted of a NS orbiting a supergiant star. Most of the new sources are 
slow and absorbed X-ray pulsars, exhibiting a large Ny{ and long Pgpin (~ 1 ks). 

• The INTEGRAL satellite has revealed the existence in our Galaxy of two previ- 
ously hidden populations of high-energy binary systems. First, the SFXTs, ex- 
hibiting brief and intense X-ray flares - with a peak flux of 1 Crab during 1-lOOs 
every ~ 100 days -, which can be explained by accretion through clumpy winds. 
Secondly, a previously hidden population of obscured and persistent sgHMXBs, 
composed of supergiant companion stars exhibiting a strong intrinsic absorption 
and long Pspin» with the NS deeply embedded in the dense stellar wind, forming 
a dust cocoon enshrouding the whole binary system. 

Apart from these observational facts, has INTEGRAL allowed us to better under- 
stand SgHMXBs and other populations of HMXBs? Do we better apprehend the accre- 
tion processes in HMXBs in general, and in sgHMXBs in particular, and what makes 
the fast transient flares so special, in the context of the clumpy wind model, and of 
the formation of transient accretion discs? In summary, has the increased population 
of supergiant HMXBs allowed us a better knowledge of these sources, compared to 
the ones that were already known before the launch of INTEGRAL? The answer to all 
these questions is probably "not yet", however we now have in hand more sources, 
and therefore more constraints to play with. Studying these populations will provide a 
better understanding of the formation and evolution of short-living HMXBs, and study 
accretion processes. Furthermore, stellar population models now have to take these 
objects into account, to assess a realistic number of high-energy binary systems in our 
Galaxy. 
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